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Abstract. New opacity distribution functions (ODFs) for several metal- 
licities have been computed. The main improvements upon previous 
ODFs computed by Kurucz (1990) are: (1) the replacement of the so- 
lar abundances from Anders & Grevesse (1989) with those from Grevesse 
& Sauval (1998); (2) the replacement of the TiO lines provided by Ku- 
rucz (1993) with the TiO lines from Schwenke (1998), as distributed 
by Kurucz (1999a); (3) the addition of the H2O lines from Partridge 
& Schwenke (1997), as distributed by Kurucz (1999b); (4) the addition 
of the H I-H I and H I-H + quasi-molecular absorptions near 1600 A and 
1400 A computed according to Allard et al. (1998). Other minor im- 
provements are related with some changes in a few atomic and molecular 
data. New grids of ATLAS9 model atmospheres for T e g from 3500 K 
to 50000 K and log g from 0.0 dex to 5.0 dex have been computed for 
several metallicities with the new ODFs. Preliminary comparisons of re- 
sults from the old and new models have shown differences in the energy 
distributions of stars cooler than 4500 K, in the ultraviolet energy distri- 
bution of metal-poor A-type stars, in the U— B and u— b color indices for 
T ofr <6750 K and in all the color indices for T eff <4000 K. 



1. Introduction 

Starting from 1992 Kurucz published new extended grids of model atmospheres 
for several metallicities [M/H] and several microturbulent velocities £. In each 
grid, the models cover a range in T e g and logg from 3500 K to 50000 K and 
from 0.0 dex to 5.0 dex, respectively. All the models were computed with the 
ATLAS9 code, which handles line opacity with the opacity distribution function 
method (ODFs). 

Later on, Castelli extended Kurucz's models by computing some more grids 
for metallicities with enhanced a elements abundances (+0.4 dex) and for £ 
values different from those adopted by Kurucz. Furthermore, she recomputed 
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the cool models in a few Kurucz grids by assuming the overshooting option for 
the convection switched off (Castelli, Gratton & Kurucz, 1997). 

Both Kurucz and Castelli models were computed by using Kurucz's (1990) 
ODFs. They are all based on the solar abundances from Anders & Grevesse 
(1989), except for iron which is not the same in all the ODFs. In fact, the 
iron abundance log(Np e /^tot) is —4.53 dex from Holweger et al. (1995) when 
ODFs for enhanced a elements abundances were computed and it is —4.37 dex 
from Anders & Grevesse (1989) when ODFs for no enhanced a elements were 
generated. Therefore, it is not straightforward to state the net effect of the a 
elements abundances on the model structure. 

Finally, the number of models is different in the various grids and also in 
some grids some models have 72 and some others have 64 plane parallel layers. 
These differences may give problems when interpolations in the models and in 
the grids are performed. 

2. The Aim of this Work 

The aim of this work is the computation of grids of homogenous ATLAS9 models, 
where "homogenous" means that: 

(1) The number of models is 476 in all the grids. Table 1 shows which are the 

models. 

(2) All the models have the same number of 72 plane parallel layers from 

lo 

E T Ross — —6.875 to +2.00 at steps of AlogT^ oss — 0.125. 

(3) All the models are computed with the updated solar abundances from 

Grevesse & Sauval (1998) and therefore with updated ODFs, which now 
also include H2O lines. 

(4) All the models are computed with the convection option switched on and 

with the overshooting option switched off. Mixing-length convection with 
1/H P =1.25 is assumed for all the models. The convective flux decreases 
with increasing T e g and it naturally disappears for T e g of the order of 
9000 K. 

3. The new Opacity Distribution Functions 

The line opacity is included in the new ATLAS9 models through new-ODFs 
based on the Grevesse & Sauval (1998) (GS98) solar abundances. Table 2 lists 
the differences between GS98 abundances and those from Anders & Grevesse 

(1989) which were used in the Kurucz (1990) ODFs. 

New-ODFs are computed for 57 values of temperature T (from 1995 K to 
199526 K, i.e. for logT ranging from 3.30 to 5.30 at steps of 0.02 dex) and 25 
values of gas pressure logP s (from —4.0 to +8.0 at steps of 0.5 dex), while Kurucz 

(1990) ODFs were computed for 56 values of T (from 2089 K to 199526 K, i.e. 
for logT ranging from 3.32 to 5.30 at steps of 0.02 dex) and 21 values of logP s 
(from —2.0 to +8.0 at steps of 0.5 dex). The larger ODF table should allow 
better modelling the upper layers of cool stars and of giant stars. 
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Table 1. The 476 models of each grid. 
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Table 1. cont. 
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Except for TiO, new-ODFs are computed with the atomic and molecular 
line lists as in the Kurucz (1990) ODFs, but some molecular lists have been 
extended to more bands (i.e. CN, OH, and SiO) and some bugs in the line lists 
have been corrected. 

In new-ODFs, the TiO lines computed by Kurucz (1993) with semi-empirical 
methods have been replaced by the TiO lines from Schwenke (1998), as dis- 
tributed by Kurucz (1999a) on the CD-ROM No 24. Furthermore, the H 2 
lines from Partridge & Schwenke (1997), as distributed by Kurucz (1999b) on 
the CD-ROM No 26, were added to the other lines. New-ODFs also include 
H I-H + and H I-H I quasi-molecular absorptions at 1400 A and 1600 A, respec- 
tively, due to the collision of neutral hydrogen in the ground-state with ionized 
hydrogen and other neutral hydrogen atoms (Castelli & Kurucz, 2001). These 
absorptions are observed as broad features in A Boo stars and in metal-poor 
A-type stars. 

The treatment of the overlapping lines at the end of the term series was 
also improved in the new-ODFs. 
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Table 2. Solar abundances different in Anders & Grevesse (1989) 
(AG89) and in Grevesse & Sauval (1998) (GS98). The abundances are 
given as log(N e ; em /H g ). 
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Figure 1. Comparison of the ultraviolet flux from the old-ODF model 
(thin dashed line) and the new-ODF model (thick full line) with pa- 
rameters T eff =8000 K, logg=2.5, [M/H]=-2.5, £=2.0 km s" 1 



New-ODFs are computed on an Alpha Digital Personal Workstation run- 
ning OPEN VMS using the codes XNFDF and DFSYNTHE. The XNFDF code 
pretabulates number densities and continua opacities. The DFSYNTHE code 
computes a line opacity spectrum for each T,P gas couple from 89.7 A to 100000 A 
at 500000 resolution. Line opacities are then sorted in the same wavelength in- 
tervals as in the Kurucz (1990) ODFs. 



4. The New Grids and Future Work 

New full grids of model atmospheres computed up to now are those for: 
[M/H]=0.0, -0.5, -0.5a, -1.0, -1.5 and £=2 km s" 1 

where "a" means that the abundances of the a elements O, Ne, Mg, Si, S, Ar, 
Ca, and Ti are enhanced by +0.4 dex over the solar or solar-scaled abundances. 
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Figure 2. Comparison of fluxes from the old-ODF model (thin 
dashed line) and the new-ODF model (thick full line) with parame- 
ters T efr =4250 K, logg=1.5, [M/H]=-0.5, £=2.0 km s _1 . The bands 
of CN and CO are indicated in the figure. 
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Figure 3. Comparison of fluxes from the old-ODF model (thin 
dashed line) and the new-ODF model (thick full line) with parame- 
ters T eff =3500 K, logg=5.0, [M/H]=0.0, £=2.0 km s _1 . The bands of 
TiO, CN, H2O, and CO are indicated in the figure. 
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Figure 4. Comparison of fluxes from the old-ODF model (thin 
dashed line) and the new-ODF model (thick full line) with parame- 
ters T efr =3500 K, logg=0.0, [M/H]=0.0, £=2.0 km s _1 . The bands of 
TiO, CN, OH, and CO are indicated in the figure. 




400 440 480 520 560 600 640 680 720 760 800 840 880 920 960 1000 1040 

Figure 5. The same as Fig. 4, but for the smaller wavelength region 
400-1120 nm plotted on a larger wavelength scale. 
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Figure 6. T e g-color relations for logg=4.5 and [M/H]=0.0 when old- 
ODF models (thick line) and new-ODF models (thin line) are used: (a) 
for U-B; (b) for u-b; (c) for B-V; (d) for b-y. 



Together with the models we have computed corresponding grids of fluxes 
and color indices in the UBV and RIJKL Johnson system, in the VRI Cousins 
system, and in the uvby Stromgren system. The new grids are available at 
http: / /kurucz. harvard.edu/grids/gridxxxodfnew. 

Small new-ODF grids for a few values of T e g and logg have been also com- 
puted for [M/H]= -1.5a, -2.0, -2.0a, -2.5, -2.5a, and -3.0a (Castelli & 
Cacciari, 2001). We plan to extend these small grids to all the models listed in 
Table 1. In addition to models and fluxes computed for £=2 km s _1 we plan to 
compute models and fluxes for £=0, 1, 4, and 8 km s _1 . Color indices in several 
photometric systems will be also computed. 



5. Comparison between Old and New Models 

Preliminary comparisons of energy distributions have shown that the different 
solar abundances do not affect the models very much. In fact, the differences 
between old-ODF and new-ODF models are significant only when additional 
absorptions introduced in the new-ODFs play an important role. 

Remarkable differences which increase with decreasing gravity and/or metal- 
licity can be observed in the 1250-2000 A region for T e g- between 9000 K and 
7000 K. They are mostly due to the H I-H + and H I-H I quasi-molecular ab- 
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sorptions. Fig. 1 shows the energy distributions from old-ODF and new-ODF 
models for T eff =8000 K, logg=2.5, [M/H]=-2.5, £=2 km s _1 . 

The new TiO line list, the addition of H2O lines and some improvements 
in the molecular data have greatly modified the energy distributions of models 
cooler than 4500 K. Examples are given in Figs. 2, 3, 4, and 5. Figure 5 shows 
in more detail the TiO and CN absorptions in the range 4000-11200 A for the 
model T efr =3500 K, logg=0.0, [M/H]=0.0, £=2 km s _1 . 

Color indices U-B and u-b from new-ODF models are redder than the indices 
from the old-ODF models for T e g<6500 K (Fig. 6a, b). The difference increases 
with decreasing T e g. For instance, the difference for U— B amounts to 75 K 
when the model T e g=5750 K, logg=4.5, [M/H]=0.0, £=2 km s~ x is considered. 

Both (B-V) and (b-y) from new-ODFs are redder than those from old-ODFs 
for T c ff<4500 K, but the differences for (b-y) are almost negligible (Fig. 6c,d). 
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